Electrodynamic Casimir Effect in a Medium-Filled Wedge II 
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We consider the Casimir energy in a geometry of an infinite magnetodielectric wedge closed by a 
circularly cylindrical, perfectly reflecting arc embedded in another magnetodielectric medium, under 
the condition that the speed of light be the same in both media. An expression for the Casimir 
energy corresponding to the arc is obtained and it is found that in the limit where the reflectivity of 
the wedge boundaries tends to unity the finite part of the Casimir energy of a perfectly conducting 
wedge-shaped sheet closed by a circular cylinder is regained. The energy of the latter geometry 
possesses divergences due to the presence of sharp corners. We argue how this is a pathology of 
the assumption of ideal conductor boundaries, and that no analogous term enters in the present 
geometry. 

PACS numbers: 42.50.Pq, 42.50.Lc, ll.lO.Gh 



The Casimir effect [l[ may be understood as an effect 
of the fluctuations of the quantum vacuum. Casimir's 
original geometry involved two infinite and parallel ideal 
metal planes which were found to attract each other with 
a negative pressure scaling quartically with the inverse 
interplate separation. In a seminal paper, Lifshitz gen- 
eralised Casimir's original calculation to imperfectly re- 
flecting plates ^ . Since its feeble beginnings research on 
the Casimir effect has grown from being of peripheral in- 
terest to a few theorists to a bustling field of research both 
experimental and theoretical with publications number- 
ingin the hundreds each year. Recent reviews include 
11, ii. 

Progress on Casimir force calculations for other geome- 
tries has been slower in coming. Spherical and cylindrical 
geometries have naturally been objects of focus, the lat- 
ter of direct interest to the effort reported herein. Only 
in 1981 was the Casimir energy of an infinitely long per- 
fectly conducting cylindrical shell calculated [Tj and the 
more physical but also significantly more involved case of 
a dielectric cylinder was considered only in recent years 
B im El El El, El- We might also mention re- 
cent work on the cylinder defined by a i5-function poten- 
tial, a so-called semitransparent cylinder [l^; for weak- 
coupling, both the semitransparent cylinder and the di- 
electric cylinder have vanishing Casimir energy. 

Closely related to the cylindrical geometry is the in- 
finite wedge. The problem was first approached in the 
late seventies [l^, E3| as part of the still ongoing debate 
about how to interpret various divergences in quantum 
field theory with sharp boundaries. Since, various em- 
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bodiments of the wedg e have been treated by Brevik and 
co-workers flsL Tgl . l2d ] and others [U, |23] ■ A review may 
be found in [23] ■ A wedge intercut by a cylindrical shell 
was considered by Nesterenko and co-workers, first for a 
semi-cylinder [24| . then for arbitrary opening angle [25l ]. 
and the corresponding local stresses were studied by Sa- 
harian [H,!!^,!!!]. The group at Los Alamos studied the 
interaction of an atom with a wedge |29l, previously 
investigated by Barton [31] and others [32, 33], the ge- 
ometry realised in an experiment by Sukenik et al. some 
years ago [s^l • A recent calculation of the Casimir energy 
of a magnetodielectric cylinder intercut by a perfectly re- 
fiecting wedge filled with magnetodielectric material was 
recently reported by the current authors [ssl ]. Common 
to all of these theoretical efforts is the assumption that 
the wedge be bounded by perfectly conducting walls. 

While until recently relatively few treatments of the 
vacuum energy of the wedge existed, the problem of cal- 
culating the diffraction of electromagnetic fields by a di- 
electric wedge within classical eletromagnetics is an old 
one and several powerful methods have been developed 
within this field. The Green's function of the potential 
(Poisson) equation in the vicinity of a perfectly conduct- 
ing wedge was found more than a century ago by Macdon- 
ald [36| and extended to the wave equation with a plane 
wave source by Sommerfcld 37]. Generalising Sommer- 
feld's method, the first theoretical solution to the scatter- 
ing problem involving a wedge of finite conductivity was 
found by Malyuzhincts in his PhD work (see [s^] for 
a review; cf. also [40|). 

A different method was proposed by Kontorovich and 
Lebedev in 1938 [ilj and used by Oberhettinger to solve 
the Green's function problem some time later The 
method has been given attention in recent analytical and 
numerical studies of the diffraction problem [il, IH, |45| . 

MM- 

In the present effort we study the Casimir energy in 
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a magnetodielectric wedge of opening angle a inside and 
outside a perfectly conducting cylindrical shell of radius 
a — See Fig. [1] The interior and exterior of the wedge 
are both filled with magnetodielectric material under the 
restriction of isorefractivity (or diaphanousness), that is, 
the index of refraction n^{uj) = e{Lu)^(uj) is the same ev- 
erywhere for a given frequency. This condition is adopted 
because without it the problem is no longer separable and 
not readily solvable. Moreover, we suspect that nondi- 
aphanous media will lead to divergences, at least in the 
absence of dispersion. 

As a natural extension of the considerations in [s^ we 
derive an expression for the free e nerg y of such a system 
by use of the argument principle [48|. (By free energy, 
we mean that bulk terms not referring to the circular 
arc boundary are subtracted.) The necessary dispersion 
relation provided by the electromagnetic boundary con- 
ditions at the wedge sides is derived in two different ways; 
by a standard route of expansion of the solutions in Bessel 
function partial waves and by use of the Kantorovich- 
Lebedev (KL) transform. (Still a third method, based on 
the Green's function formulation, is sketched in the Ap- 
pendix.) The corresponding boundary condition equa- 
tion at the cylindrical shell is well known. These together 
allow us to sum the energy of the eigenmodes of the ge- 
ometry satisfying eigenvalue equations for the frequency 
and azimuthal wave number v by means of the argument 
principle. 

There are important differences between the di- 
aphanous geometry considered herein and the standard 
geometry of a perfectly conducting wedge. Assuming di- 
aphanous electromagnetic boundary conditions, the in- 
terior and exterior wedge sectors are coupled and re- 
main so also in the limit where the reflectivity of the 
wedge boundaries tends to unity (for example by letting 
e — > oo, /i — > so that their product is constant). Assum- 
ing the wedge be perfectly conducting from the outset, 
however, the interior of the wedge is severed cleanly from 
its exterior at all frequencies, a significantly different sit- 
uation. 

The Casimir energy of the perfectly conducting wedge 
and magnetodielectric arc considered in 35] was found to 
possess an unremovable divergent term associated with 
the corners where the arc meets the wedge. This is a 
typical artifact of quantum field theory with non-flat 
boundary conditions (e.g. [2l|, HH])- We will argue 
in section IIIBI that there is no such term present in 
the geometry considered herein, and that the direct gen- 
eralisation of the finite part of the energy of the sys- 
tem considered in |35|] to the present system is in fact 
the full regularised Casimir energy. The reason for this 
rests upon two unphysical effects of perfectly conduct- 
ing boundary conditions at the wedge sides (the van- 
ishing of the tangential components of the electric field 
there). Namely, such boundary conditions exclude the 
existence of an azimuthally constant TM mode, and di- 
vides space cleanly into an interior and exterior sec- 
tor with no coupling allowed between modes in the two 
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FIG. 1: The wedge geometry considered. 



sectors. Moreover, for a wedge consisting of perfectly 
conducting thin sheets dividing space into two comple- 
mentary wedges, the ideal conductor boundary condi- 
tions will count the azimuthally constant TE mode twice 
whereas with more realistic boundary conditions such as 
considered here, such a mode must be common to the 
both sectors, < 9 < 2tt. In these two respects the 
perfectly conducting wedge differs from the diaphanous 
one, and put together these redefinitions provided by the 
diaphanous wedge exactly remove the divergent extra en- 
ergy term found in [35] and previously in |253 . 

We show numerically that except for the singular term, 
the energy of a perfectly conducting wedge closed by 
a magnetodielectric cylinder whose reflectivity tends to 
unity is regained in the limit where we let the wedge 
boundaries become perfectly reflecting. 

I. BOUNDARY CONDITIONS AND 
DISPERSION RELATIONS 

We begin by considering in general the form of an ex- 
pression of the energy of a diaphanous wedge inside and 
outside a cylindrical shell such as depicted in Fig. [1] We 
assume the cylindrical shell to be perfectly reflecting. Let 
the interior sector —a/2 < 9 < a/2 have permittivity and 
permeability ei and /ii relative to vacuum, and the cor- 
responding values for the exterior sector tt > \9\ > a/2 
be €2 and fi2 so that ei(a')/zi(w) = €2(10) ^2{^^) = n^{(^)- 
The cusp of the wedge is chosen to lie along the 2 axis, 
which is also the center of the cylindrical shell, and the 
interfaces are found a,t 9 = ±a/2 and at p = a (p is the 
distance to the z axis). 

We will calculate the Casimir energy by 'summing' over 
the eigenmodes of the geometry using the so-called argu- 
ment principle, now a standard method in the Casimir 
literature. The eigenmodes of a given geometry are given 
by the solutions of the homogeneous Helmholtz equation 

{V"^ -n^df)u{r,t) ^0, (LI) 

which also satisfy the system's boundary conditions. 
Here u symbolises a chosen field component of either 
the electric or magnetic field. We will choose and 
Hz as the two independent field components from which 
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the rest of the components can be derived by means of 
Maxwell's equations. 

The translational symmetry with respect to z and time 
makes it natural to introduce the Fourier transform 



where p = (/o, ff) and p = \J + y^- The Helmholtz 
equation now simphfies to the scalar Bessel equation 

(Vi + fci)i?,(p;fc„c^) = 0, (1.2) 

where 

^\=dl + -d, + \dl (1.3) 
P P 

and k\ — e^w^ - fc^. 

We will define the quantity k as 

K = ■\/fc2 — e/iw^ — —ikj_, (1-4) 

where the root of k is to be taken in the fourth com- 
plex quadrant. When in the end we take frequencies to 
lie on the positive imaginary axis, k becomes real and 
positive, something we bear in mind in the subsequent 
calculations. 

A general solution to Eq. p.2p is of the form 

- [AM'\k^p) + B,Mk^p)]iaeJ'''+be~''"'), (1.5) 

where A^, B^, a, b and v are arbitrary. If, as in our case, 
ly is allowed to take non-integer values, we must restrict 
it to > because except at integers J^{z) and J-^{z) 
are linearly independent. 

Solutions of the EM field in a wedge geometry are ex- 
pressed as a sum over cylindrical partial waves whose 
kernel are Bessel and Hankel functions of argument k^p. 
Thus it is clear that the boundary conditions on the 
wedge surfaces can only be solved for each partial wave 
if the speed of light is the same in both sectors since fcj^ 
would otherwise take different values in the two media for 
given kz and oj and the kernel functions would be linearly 
independent functions of these. The diaphanous condi- 
tion is thus prerequisite for the explicit solution of bound- 
ary conditions below. Without this condition the prob- 
lem at hand is not analytically solvable with the methods 
used herein. We expect that even if we could solve the 
nondiaphanous problem we would encounter divergences 
that might or might not be curable by the inclusion of 
dispersion. 

The presence of the wedge primarily has the role of 
dictating which values of v are allowed. If one were to 
consider a cylinder (periodic boundary conditions), only 
integer values of v (both positive and negative) would 
be acceptable and expressing the solution as a sum over 
these integer values would be appropriate. Were one in- 
stead to let the wedge be perfectly reflecting (Dirichlet 
and Neumann boundaries at ±a/2 where a is arbitrary) 



V would be forced to take values that are non-negative in- 
teger multiples of 7r/a. The diaphanous magnetodielec- 
tric boundaries present here also restrict v to discrete 
values for given e's and /i's, but explicitly determining 
these values is no longer immediate because modes ex- 
isting in the exterior and interior sectors now couple to 
each other. For a given frequency we therefore make use 
of an appropriate dispersion function representing these 
boundaries in order to sum over the appropriate values 
of V by means of the argument principle, whereupon we 
may sum over the eigenfrequencies of the modes inside 
and outside the cylindrical shell to obtain the energy. 

The boundary condition dispersion relation pertaining 
to the circular boundary is known (e.g. [s^, Eq. (4.12), 
with = 1), 

g,{k,,^) = l-x^\l{x)=Q, (1.6) 
where x = an, 

K{x) = ^[I,{x)K,{x)l (1.7) 
ax 

1 1, , Ki, are the modified Bessel functions of the first and 
second kind of order v. We can simply use this equation 
to sum modes satisfying the boundary condition on both 
sides because the wedge boundaries at iba/2 impose the 
same discretisation of v inside and outside the cylindrical 
shell (were we to have e.g. a third medium in the sector 
1^1 < Oi/2,p > a different from medium 1, this would 
no longer be the case as we will see: the field solutions 
would take different values of v inside and outside the 
cylindrical boundary and the boundary conditions at the 
cylinder could no longer be solved for each eigenvalue of 
v) . We now turn to a derivation of the dispersion relation 
pertaining to the interfaces at = ±q;/2. 

In the following we shall use the term TE to denote 
electromagnetic modes whose E field has no component 
in the z direction and TM denotes those modes whose H 
field has no z component. This is not 'transverse elec- 
tric' and 'transverse magnetic' with respect to the wedge 
boundaries at 6* = ±q;/2, but this does not matter since 
we will find that the eigenequation of these boundaries 
is the same for all field components by virtue of the di- 
aphanous condition. 

A. Kontorovich-Lebedev approach 

We will first employ the technique of the Kontorovich- 
Lebedev (K-L) transformation [4l| and its inverse trans- 
form which may be written as: 

/>ZOO 

Ez{p) —i I dviye^ sm{Trh')Ki,{Kp)<iz{9;h'), (1.8a) 

^0 

t,{9-v)^— -^e~'^K,{Kp)E,{p), (1.8b) 
Jo P 

(dependence on k^ and uj is implicit). While less exten- 
sively covered in the literature than most other integral 
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transforms, some tables of K-L transforms exist f49|, |50|. 
Numerical methods for evaluating such transforms were 
recently developed by Gautschi [5l|. We will ignore the 
presence of the cylindrical shell in this section and only 
study how the presence of the walls of the wedge dis- 
cretize the spectrum of allowed values of the Bessel func- 
tion order v. 

With this, (|1.2p . after multiplying with p^, transforms 

to 



(1.9) 



Equation (|1.9p is now in a form fully analogous to that 
encountered in a planar geometry fe. g. Is^ . [ssl. [53|V We 
follow now roughly the scheme of I54| and determine 
the dispersion relation [condition for eigensolutions of 
Eq. (|1.2p ] by means of summation over multiple reflec- 
tion paths. By noting that the solutions to Eq. (|1.9p have 
the form of propagating plane waves travelling clockwise 
or anticlockwise along the now formally straight 9 axis 
{u playing the role of a reciprocal azimuthal angle) the 
analogy to a plane parallel system is obvious. 

We write the solution of Eq. (|1.9p in the interior sector, 
\9\ < a/2, in the form 



e e 



-fee' 



(1.10) 



where e^ are undetermined integration coefficients, field 
amplitudes at = to be determined from boundary 
conditions at — ±a/2. 

Likewise the solutions in the exterior sector (the 'com- 
plementary wedge') TT > 16*1 > a/2 may be written 



(1.11) 



where the undetermined amplitudes e^ are 'measured' at 
= 11. The choice to measure the amplitudes in sectors 1 
and 2 at = and tt respectively is arbitrary, but makes 
for maximally symmetric boundary equations. 

The homogeneous Helmholtz equation thus solves the 
scattered part of <tz given some source field <!f°. Let us 
assume there is a source field in form of an infinitely thin 
phased line source parallel to the z axis at some position 
6*0 in the interior sector. The direct field (which only 
propagates away from the source) may be written in the 
form 



6(61 



0(00 - e)e-^ 



-iv9 



(1.12) 



where Q{x) is the unit step function and the field ampli- 
tudes are "measured" aX 9 = We do not need to know 
the constants ej explicitly and take these to be known 
constants. The multiple reflection problem (or equiva- 
lently, boundary condition problem) is now a system of 
four equations for the four amplitudes e^ , e^ as functions 
of ej. 

We define the reflection coefficients at the boundaries 
9 = ±a/2 as the ratio of reflected vs. incoming field 
amplitude, r = <tz,icfi/'^z,in, as seen by a wave coming 
from and refiected back into sector 1 (a wave going the 
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FIG. 2: The solutions of the dispersion relation (|1.20p as a 
function of r and v for a = 0.75. The eigenvalues of u for a 
given r are marked; the energy is calculated by summing over 
these values and then integrating over all uj. 



opposite way experiences a coefficient — r). With the 
assumption ei/ii = £2/^2 the reflection coefficients of the 
s and p modes differ only by a sign: 



£2 - £1 
£2 + £1 



^J■2 - Ml 



(1.13) 



We will simply use r in the following, representing either 
of the modes. We also define the transmission coefficient, 
the ratio of the transmitted to the incoming amplitude, 
going from sector i to sector j, tij where i,j = 1,2 de- 
notes the sectors in figure [TJ 



2£,, 



Mi + 



— tij.p — 



2nj 



fij + Hi 



(1.14) 



Since these coefficients are invariant under K-L transfor- 
mation, they are the sought-after single interface reflec- 
tion and transmission coefficients also in the K-L regime. 
Note that with the diaphanous condition, reflection co- 
efficients are independent of v, something which would 
not be true in general. Were r to depend on v this would 
give rise to corrections to the energy expression derived 
in section III Al (See also the Appendix, where such v 
dependence does occur.) 

We formulate the electromagnetic boundary conditions 
in terms of refiection and transmission. In the K-L do- 
main the system looks and behaves analogously to the 
planar system (see [Hll for details on this case) , but with 
one important difference, namely that a 9 directed partial 
wave which is transmitted at a wedge boundary does not 
disappear from the system but is partly transmitted back 
into sector 1 again circularly. Thus we obtain four equa- 
tions for the four amplitudes e"'", e~, e^ and e^, coupling 
to each other through paths refiected or transmitted at 
one interface: 



+ i2ie e 



0+ G^^'^ 



=ren e 



e =11260 e 



re e 
ti2e ' e" 
ti2e^e* 



(1.15a) 
(1.15b) 



s+ -U^e+e"''' + ti2e+e''"^ - rg-e'^^^^""), (1.15c) 



(1.15d) 
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Eigenvalues of v for the wedge correspond to solutions 
of these boundary conditions, which exist when the sec- 
ular equation of the set of linear equations for and 
is fulfilled. The characteristic matrix is 

/ 1 -re*''" -^216''''' \ 
_ -re''^" 1 -^216*'''' 

- -il2e''''' 1 re'«^(27r-a) 

\ -il2e*'''^ ^gii'(27r-a) ^^ J 

(1-16) 

and the dispersion relation sought after is 

= detD = 0. (1.17) 

The matrix form p.l6p is rather instructive. Note how 
D is a block matrix of the form 

Gi2 D2' ) ' 

where describes multiple scattering within sector i 
and Gij describes coupling between the sectors by trans- 
mission from sector itoj. Since the G matrices commute 
with the D matrices, det D can be written 

det D = det(DiD2 - G21G12). (1.18) 

We may use the energy conservation relation 

^12^21+^2 = 1, (1.19) 

together with (jl.lSp to find the simple expression 

D{iy,Lu) =(1 - e2™)2 - r2(g«^(2^-a) _ ^^^.a■^2 

^ - 4e2^*''[sin2(iy7r) - sm'^{i,{T: - a))]. 

(1.20) 

It is noteworthy that this dispersion relation only has 
an implicit dependence on to through the quantity r'^{uj). 
As an example we plot the solutions to Eq. (|1.17p as a 
function of v and r in Fig. [2] for a = 0.75 radians. 

Note at this point that whenever r is real, all zeros of 
D{v,uj) in (|1.20p are real. In the following we shall think 
of r as well as the eigenvalues of v as real quantities. For 
real frequencies lo reflection coefficients will in general be 
complex, while after a standard rotation of frequencies 
onto the imaginary frequency axis these coefficients are 
always real as dictated by causality. Although zeros are 
complex the argument principle may still be used; the 
discussion of connected subtleties may be found in e.g. 
0,111111]. 

It is easy to see that this dispersion relation gener- 
alises that for a cylinder (of infinite radius) and a per- 
fectly conducting wedge. In the latter limit, r = 1, the 
determinant detD has zeros where v = rrvKja and at 
V = rmr/{2'iT — a) where m is an integer. This becomes 
obvious when noting that 

D{v, uj) -Ae^'^'^" sin va ■ sin v{2ti - a). (1.21) 



This reproduces, in other words, the case where the 
wedge is made up of thin perfectly conducting sheets. 
For the perfectly conducting wedge it is customary to re- 
strict V to values that are integer multiples of 7r/a from 
the beginning. 

Likewise when the two materials become equal, 

D{v, Lo) ^ -46^'^*'' sin^ vtt = Do{u), (1.22) 

which has double zeros where v = a positive integer, 
corresponding to a clockwise and an anticlockwise mode 
or, if the reader prefers, the sum over v = +m and — m. 
This is just the cylinder case 0, i, i, [13, [nl, [13, [H, [ll- 
We see from Fig. [5] that except for i/ = which remains 
degenerate, the double zeroes split into two separate sim- 
ple zeros for finite r. For special opening angles which are 
rational multiples of tt there will be other zeroes which 
remain degenerate as well. 

One sees directly that were we to solve Eq. (|1.2[) for Hz 
instead of Ez the dispersion relation would be identical 
to Eq. (|1.17p since the only difference would be the sign 
of the reflection coefficient (we would employ Ts rather 
than Tp), which only enters squared. One should note 
that the distinction between and rp here does not cor- 
respond to the distinction between TE and TM modes of 
the entire cavity, but this is of no consequence in the fol- 
lowing because the dispersion relation (|1.17p is the same 
for all field components! 



B. Derivation by standard expansion 

We will now sketch how the result (ll.20p may be de- 
rived by a more standard method similar to that made 
use of in [35] . The solutions of Eq. (|1.5p that correspond 
to outgoing waves at p ^ cxd may be expanded following 
the scheme of [s^ in an obvious generalisation of those 
found in [57^]. Due to criteria of outgoing- wave bound- 
ary conditions at p — > 00 and non-singularity at the ori- 
gin the solution must consist purely of Hankel function 
H^\k±p) far from the origin and only of terms contain- 
ing J^{k±p) near p = 0. Following the scheme of [35,] we 
choose Hj)''\kj_p) for p > a and J,j{k±p) for p < a both 
in the interior sector —a/2 < 6 < a/2 and outside, and 
couple the solutions across these straight boundaries. It 
will not matter which Bcssel function we choose for the 
present purposes: the resulting solution expansions are 
identical but for the replacement of one Bessel function 
with another. 

In a straightforward generalisation of the expansion 
used in [ssj we write down the following general solutions 
in sector 1 of Fig. [1] for p > a: 
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Er,l — 

Ee,i = 
Ez,i — 

HrS = 

He,i = 
Hz,i ~ 



dv 



ik. 



A: I "'^ kip 



dv 



klp' 



cos vO + 1 
cos vO + i 



poo 

/ di^Hl^'> [ai cos 1^6 + ia^ sin i^9] e"" 
Jo 



dv 



dv 



k^p k±_ 



/2 

COS u6 -\- i 



/ dvHl^'> [bi cos ve + ib^ sin vO] e'"" 
Jo 



cos vO + i 

/2 



vk 



sm vU > e 



■vl2 

1 

sm e"''^^ 



K^p k 



k^ -1 fcip^" 



(1.23a) 
(1.23b) 
(1.23c) 
(1.23d) 
(1.23e) 
(1.23f) 



r 



where we have omitted the arguments of Hl^\k±p) and 
its derivative, and of di{v),bi{v), etc., the latter being 
undetermined coefficient functions of v. 

We write the solution in sector 2 in exactly the same 
form but with the simple replacements ^ — tt and 
£1 — ^ e2,Mi ~* M2 and the same for the coefficient func- 
tions. With the isorefractive assumption k± is the same 
in both media for given u; and kz, so the boundary con- 
ditions at the interfaces can be solved under the integral 
signs. In general there are 8 unknown functions and eight 
equations, yet one finds that the s and p modes decouple 
into linear equation sets of 4 x 4 on the form 



where D equals 



sm ■ 



D 







(1.24) 



-eism^ 




-ei cos ^ 



62 sini^(| - tt) 




tt) 
— sin v(j 





£2 C0Si/(f — tt) 



and a is a vector, either (ai, fl]^, 02, flj) or (61, 6]^, 62,^2 )• 
As before the eigenmodes of the system solve the equa- 
tion det D = 0. With some manipulation we find that the 
determinant can be written simply as 



det D =^(£2 — £1)^ sin^ v{Tr — a) 
- ^(£2 + £1)^ sin^ VTT. 



(1.25) 



Under the assumption that £2 + £1 7^ the equation 
detD = is equivalent to Eq. IlTTI) with Eq. (fOO]) . 



II. THE CASIMIR ENERGY 

In order to find the Casimir energy we shall employ 
the argument principle, introduced to the field of Casimir 
energy by van Kampen et al. who rederived Lifshitz's 



result in a simple way. For a very readable review of the 
technique, see [ssj . 

A similar system to that shown in Fig.[T]was considered 
in [ssj where the plane sides of the wedge were instead 
made up of perfectly refiecting interfaces and the circular 
boundary was diaphanous. We will start from the result 
of jssjl and generalise this step by step to approach the 
desired energy expression for the current situation. Ex- 
cepting the formally singular energy term associated with 
the sharp corners where the arc meets the wedge walls 
found in that paper (we shall regard this term separately 
below) , the Casimir energy per unit length of that system 
in the limit of perfectly reflecting circular arc was (Eq. 
(2.10) or (4.11) of [H]) 

(2.1) 

with g^{kz,uj) given in Eq. (|1.6[) and we deflne the short- 
hand 



TT 



(2.2) 



The prime on the summation mark means the to = 
term is taken with half weight. The integration contour 
A is chosen to follow the imaginary axis and is closed to 
the right by a large semicircle thus encircling the positive 
real axis. The roots of (|1.6p are in general complex; the 
applicability of the argument principle for such situations 
was discussed in 0, [53, HI, H^l ■ The energy S has been 
normalised so as to be zero when the circular arc is moved 
to infinity. 

Each frequency satisfying gmp{kz,uj) — gives a pole 
which adds the zero temperature energy ^ of that mode 
through Cauchy's integral theorem. In the end there are 
sums over the eigenvalues of v, mp, the eigenvalues found 
when the sides of the wedge are assumed to be perfectly 
conducting. Employing such an assumption from the 
start completely decouples the interior sector \6\ < a/2 
from the exterior. If we were to interpret the perfectly 
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reflecting wedge as the limit of an isorefractive wedge 
such as that described by the dispersion relation (|1.20p 
as \r\ — > 1, however (for example by letting 62 — > oo 
and ^2 — > so that their product is constant), the inte- 
rior and exterior sectors remain coupled and wc obtain 
an additional m-sum, namely that over i' = mp' of the 
complementary wedge, where 



P 



27r - a 2p-l 



(2.3) 



To obtain direct correspondence we therefore modify 
Eq. (|2.ip by also including the energy of the modes of the 
complementary wedge, fulfilling v = mp' . Since we will 
soon generalise this result to the case where the wedge 
is diaphanous, it is reasonable to subtract the energy 
corresponding to absence of the boundaries at ±a/2, by 
subtracting off the energy obtained were v to fulfill peri- 
odic boundary conditions (i.e. a circular cylinder). The 
result is 



id 



1 

47ri 



dujuj ^ In "^^^ ' -^^P' 



duj 



5,1 



(2.4) 

The periodic function gm{kz,uj) is squared since both 
positive and negative integer orders contribute equally 
in the periodic case, and the symmetry under m — > —m 
makes for a factor of 2 except for m = 0. The latter 
exception is automatically accounted for by the prime on 
the sum. 

Note that employing g^{kz,uj) with the argument prin- 
ciple automatically takes care of the sum over the two 
polarisations since, by virtue of the diaphanous condi- 
tion, (7^ is a product of boundary conditions for TE and 
TM modes (see e.g. Appendix B of [sslp. 

Let us now perform the generalisation of Eq. (|2.4p to 
the present case. The sum over v = mp and mp' may be 
generalised to a sum over the solutions of Eq. (|1.17p using 
the argument principle once more to count the zeros of 
Eq. (|1.20p . and the subtraction of the periodic modes in 
the absence of the boundary is performed by subtracting 
the solutions of Dq{v) = with Dq from Eq. (|1.22p (note 
that the zeros of Dq are double, automatically giving the 
factor 2 manually introduced in Eq. (|2.4p by taking the 
square of gm)- We obtain 



1 



2(27ri)2 



dk. 



—00 
d 

duj 



dhJUJ 



2vr JA 
h\gy{kz,uj) 



Un^. (2.5) 



The contour of the v integral is the same as that for the 
Lo integral. 

Neither of the contour integrals obtain contributions 
from the semicircular contour arcs so we are left with 
integrals over imaginary order and frequency. Performing 



substitutions uj = iC, and v = i-q we obtain 
1 



S =- 



dk. 



dCC 



drj 



— In gi,,{kz,iC) 



Am^i!^. (2.6) 

dr] Do{irj) 



This is the general form of the Casimir energy of the 
system presented herein. 

To be very explicit about the regularisations per- 
formed: Eq. (12. 6p is the energy of the geometry of Fig. [U 
minus the energy when the cylinder is pushed to infinity 
(the double wedge alone), minus the renormalised energy 
of the cylinder relative to uniform space. 



S — {S© 



(2.7) 



where (So is the C-renormalised energy of a cylindrical 
shell (relative to uniform space) considered in pi)], and 
@, < symbolise the double wedge with and without the 
cylindrical shell. It is thus clear that the energy should 
vanish when either the cylindrical boundary tends to in- 
finity (^g, — > (o'< and So — > 0) or when the wedge becomes 
completely transparent (<^g, — So). 

The corresponding free energy at finite temperature 
T is found by simply substituting the integral over C in 
Eq. (|2.6p with the well known Matsubara sum over the 
frequencies Qk — 'inkT where k G Z; 



/oo 
-CXD 7 



(2.8) 



We will not consider finite temperature numerically in 
the present effort. 



A. Non-dispersive approximation 

In order to proceed to producing numerical results we 
make the simplifying assumption that r be approximately 
constant with respect to C over the important range of C 
values; ^ w 0. This a version of the constant reflection 
coefficient model which was previously found to be useful 
for the planar geometry jGO] . While it is true that for any 
real material, reflectivity must tend to zero at infinite fre- 
quency, the non-dispersive approximation is a useful one 
and allows a simpler expression to be derived. We will 
find below that the resulting Casimir energy expression 
is finite even when r — 1 for all frequencies, except when 
a = or 27r. There is therefore no need to assume high- 
frequency transparency for the sake of finitcness in this 
case. 

With this assumption we can easily perform a partial 
integration in We note that, when r is independent 
of 1/ as in the diaphanous case (see the Appendix for a 
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situation where this is not so), 



so that X'^ + Y'^ = and 



d 
drj 



In- 



D 



Q;sinhr/(27r — a) — {2tt — a) sinhrya 



sinh^ rjTT — sinh^ ri{'K — a) 
sinh77(7r — a) 
sinh r/TT ' 



(2.9) 



which is now approximated as independent of C and k^. 
It is then opportune to perform a change of integration 
into the polar coordinates 



X — = KCOs9; Y = kz ~ nsniO, 



(2.10) 



dkz / dUian) = — / dxxf{x), (2.11) 



where x = an as before. We obtain after integrating by 
parts 



drj- 



' sinh7](7r — a)[a sinh 7^(271 — a) — (27r — a) sinhyya] 
sinh ?77r [sinh^ rjiT — sinh^ r;(7r — a)] 

I 



dxx\i\\V~x \^(x)\. 



(2.12) 



Despite appearances this expression is in fact real. This 
is because the dispersion function in the first integral is 
an odd function of t] while the real and imaginary parts 
of the logarithm are even and odd respectively (provided 
the appropriate branch of the logarithm is taken) , hence 
the imaginary part of S vanishes under symmetrical in- 
tegration. It is straightforward to write down the cor- 
rection terms containing 4? or should the reader wish 
to do so. Such is necessary were one to study the role 
of dispersion on the energy; we shall not consider this 
herein — but see the Appendix for dr/dv ^ 0. 



The energy expression (|2.12p has the reasonable prop- 
erties of being zero at a = tt and symmetrical under the 
substitution a <-> 27r— a. We will study Eq. (|2.12p numer- 
ically in section lTTTl We argue in the next subsection that 
Eq. (|2.6p is the full Casimir energy of this system (after 
subtracting that of the cylinder alone). Thus the zero 
energy at a = tt demonstrates a particular generalisation 
of the theorem of Ambj0rn and Wolfram ([61|, stated in 
Eq. (2.49) of Q): the energy of a semi-circular compact 
diaphanous cylinder is half that of a full cylinder (there 
is an equal contribution from the exterior 'half- cylinder' 
so the difference is zero). 

For large rj the term proportional to a in the big frac- 
tion in Eq. (|2.12p behaves for tt — a > as 



d , D 2ar^ 
— In — 

d-q Do e^"") - 



(2.13) 



with a similar behavior for the term proportional to 
27r — a, and so is exponentially convergent. With per- 
fect reflectivity Eq. (|2.12p is finite except when a equals 
or 2tt when Irl = 1. 



B. No additional corner term 

In the geometry considered in [s^ , which differed from 
the present one primarily by the assumption that the 
wedge be perfectly conducting, the Casimir energy was 
found to possess a divergent term which could be asso- 
ciated with the corners where the arc meets the wedge 
sides. When the arc was instead made diaphanous it was 
shown that this term could be rendered finite by virtue 
of high frequency transparency as displayed by any real 
material boundary. 

The energy (|2.6p is the direct generalisation of the fi- 
nite part of the energy of the system considered in [3^ . 
We will argue that when the wedge is also diaphanous, 
this is indeed the full energy of the system, regularised 
by the subtraction of the energy of the cylinder alone 
(which in turn is regularised by subtracting the energy 
of uniform space). 

Let us recapitulate how the divergent term in [35| came 
about. The zeta function regularised energy expression 
(Eq. (4.13) of [IHl) adds the m = modes of both polar- 
isations with half weight. There should be no m = 
TM mode, however, because the perfectly conducting 
wedge forces any azimuthally constant electric field to 
have zero amplitude everywhere, thus the half-weight 
zero TM mode should be subtracted. Moreover, since 
for arbitrary opening angles only positive values of ni 
are allowed, the zero TE mode should be counted with 
full rather than half weight, and thus the correction term 
equals one half the to = energy of the TE mode minus 
one half that of the TM mode. 

In contrast we are here not considering perfectly con- 
ducting wedge boundaries so the TM to = mode should 
be included. The question becomes whether the v = 
TE and TM modes have been counted with only half the 
weight they should. In a system such as ours the interior 
and exterior sectors are coupled and all allowed modes 
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are modes satisfying boundary conditions of the whole 
double wedge. Thus there can be only one azimuthally 
constant mode for all (not one for each sector as one 
obtains for a perfectly conducting wedge-sheet) hence the 
zero mode should be counted once. This is exactly what 
is done in Eq. (12. 6p because the dispersion function (|1.20|) 
has a double zero at = cancelling the factor 1/2. 
Hence no additional correction term is necessary and the 
use of dispersion relations with the argument principle 
automatically gives the full result. 

In our numerical considerations reported in section Hill 
we find correspondence with the finite part of the en- 
ergy reported in [sHi when applied to two complemen- 
tary wedges separated by a perfectly conducting sheet. 
Note how this correspondance is somewhat peculiar: In 
the energy expression of that reference the zero mode 
was counted with half weight where it should have been 
accounted for fully, but in adding the energy of the com- 
plementary wedge as in Eq. p.9p each of the complemen- 
tary wedges contribute a half of the m = mode energy, 
amounting to the full energy when we insist that this 
mode be common to the whole system. 

It is thus made clear how the divergent term found 
in [U, m, [IHl can be seen as a pathology of the ideal 
conductor boundary conditions a.t 9 = ±a/2 which (a) 
completely removes the azimuthally constant TM mode 
and (b) cleanly severs the connection between the inte- 
rior and exterior of the wedge. Whether a similar term 
would appear - perhaps with a finite value ~ for a non- 
diaphanous wedge remains an open question, since the 
diaphanous condition employed herein is also a special 
case. 



III. NUMERICAL INVESTIGATION 

It is useful to introduce the shorthand notation 
1 



100 



47r2r 



drj 







^ In ^ 
dr/ Do 



(3.1) 



where Y is the imaginary part of the integral over the 
logarithm in Eq. ([^l^ . 



y{ri) 



dxa;arctan- 
1 -a;^5Re{A4(x)} 



(3.2) 



where we take the argument of the logarithm to lie in 

r 7r 7r 1 

Near a; = this integrand behaves like a::sin(lna;), os- 
cillating increasingly fast. Techniques of rotating the in- 
tegration path are restricted by the scarcity of methods 
for evaluating Bessel functions of general complex order, 
and will anyway come at the cost of making ^ In Z?/ Dq 
oscillatory. For numerical purposes it is more useful to 
perform the substitution x = e^: 



^iv) = / dye^y arctan 



-e2v3m{Af,(e^)} 



(3.3) 
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FIG. 3: (a)Different methods of calculation used in different 
areas of the x,r] plane (see text), (b) The function Ylrf). 



For moderate values of rj this integrand is numerically 
manageable (there are ^{Ari) significant oscillations to 
integrate over), the remaining challenge being the evalu- 
ation of Xiri{x). 

Rather than consider the complex function lir^lx) it is 
numerically useful to consider the real function 



(3.4) 



When rj is real, Lir,{x) = Ke/i,,(x). We find, using the 
Wronskian relation 

W[K^,h]{x) ^l/x (3.5) 

and relations between the two modified Bessel functions, 
that Air, can be written as 

1 / / X / \ 2i sinh rjir . 
Xiriix) = - + 2Ki Jx)Lir,[x) (x)Kin(x) 

X ' TT ' 

(3.6) 

For obtaining the right limit of the integrand near rj = 
one may notice that Y{ri) for small r/ is 



'Tj I dxx^ 



AKoKiil - 2xIoKi) 



,0 X- - (1 ~ 2x/oifl)2 

plus higher orders. Numerically one finds 

y(ry) ~ 0.87442?7 + ff{r]^). (3.7) 
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A complete algorithm for evaluating K and L for imag- 
inary order and real argument was developed by Temme, 
Gil and Segura [13, [63]. Since we are only calculating 
products of Bessel functions and the methods for calcu- 
lating one is much like that for another, the code perfor- 
mance could be increased significantly by reprogramming 
(we used programming language C#). 

Different calculation methods are appropriate in differ- 
ent areas of the x, 77 plane as shown in Fig. [3^. For K and 
K' we use Maclaurin type series expansion in region I in 
the figure (bounded by 77 > 0.044(x — 3.1)^-^-|-x — 3.1) and 
in regions 11 and 111 (bounded by 77 < 380(g^)°-^'''^) a 
method of continued fractions is used [6J] (the continued 
fractions method of [65l | may be used for imaginary or- 
ders also). No continued fractions method is available 
for L, but series expansions turn out to be more ro- 
bust than for K'; for x < 60 (region 11) Maclaurin 
series expansion is used, and asymptotic series expansion 
is used above this (region 111) . In the remaining area (re- 
gion I V) Airy function type asymptotic expansions were 
used [H, [661 \6T\ . In addition a fast method for evalu- 
ating complex gamma functions was necessary - we used 
that of Spouge [g^ . The resulting algorithm was able to 
calculate Ai^(x) with at least eight significant digits on 
x, 77 e [0, 100], more than sufficient for our purposes. 

Because the calculation of A is rather elaborate we 
do not do the double integral (|3.ip directly, but calcu- 
late a number of discrete values of Y{ri) and use spline 
interpolation to represent Y in the integration over 77, 
which then converges rapidly. The function Y{ri) is zero 
at 7/ = and increases smoothly thence to approach a 
positive constant, obtained already at modest values of 
77, as plotted in Fig. [3b- The factor [In D/Dq]' behaves 
as e"^''" for large 77 (assuming a < it) assuring rapid 
convergence when a is not close to zero or 27r. 

In the limit ^ 1 we should obtain correspondence 
with [ssj where the energy of a perfectly reflecting wedge 
closed by a diaphanous arc was considered. In this strong 
coupling case (the arc becoming perfectly reflecting) the 
energy of the sector inside the wedge only (modulo a 
singular term) was written on the form 

1 



3(P), 



(3. 



where the dimensionless function e{p) is given in [35| . Eq. 
(4.22), and p — ir/a as before. In the present case the 
modes in the interior and exterior sectors never decou- 
ple, even in the limit r ^ 1 and Eq. (|2.12p thus calculates 
the energy of the whole system, regularised by subtract- 
ing the energy of free space, that is, by subtracting the 
result when the arc is moved to inflnity (this is already 
implicitly subtracted by use of Eq. p.6p ) and the wedge 
boundaries become transparent. The energy to compare 
with is therefore on the form given in (ssj where the en- 
ergy of the complementary wedge is added and that of a 
cylinder subtracted. We can therefore write Eq. (|2.4p in 
the form (fid = eid(p)/87r7ia^ where 
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FIG. 4: Above: The function e(p) calculated for different r. 
Below: Same quantity, now plotted as a function of opening 
angle a. 



and p' was deflned in 

For our system the corresponding function is 



e{p) 



2 



C?77 



In- 



D 



Y{iri; r). 



(3.10) 



We plot e{p) as a function of p and as a function of a 
in Fig. 31 When r ^ 1 full agreement with Cidip) of 
Eq. (pJ)) is obtained. 



eid(p) = e(p)+e(p')-2e(l). 



(3.9) 



Conclusions 

We have analysed the Casimir energy of a magnetodi- 
electric cylinder whose cross section is a wedge closed by 
a circular arc under the restriction that the cylinder be 
diaphanous, i.e. that the speed of light be spatially uni- 
form. We obtain an expression for the Casimir energy per 
unit length of the cylinder, regularised by subtraction of 
the energy of the wedge alone and the cylinder alone. 
The energy is then zero when the opening angle of the 
wedge, a, equals tt, it is symmetrical under the substitu- 
tion a «-> 27r — a, and it remains finite as a tends to zero 
or 2tt except when the absolute reflection coefficients of 
the wedge boundaries equal unity. 

A numerical investigation confirms that this gener- 
alises a previously known result for a perfectly conduct- 
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ing wedge closed by a diaphanous magnetodielectric arc 
in the hmit where the arc becomes perfectly reflecting, 
except for a singular term present in that geometry which 
we argue does not present itself in the present configura- 
tion. This implies that the singular term found and dis- 
cussed in 35] is an artefact of the use of ideal conductor 
boundary conditions and does not enter for a diaphanous 
wedge. 

We mention finally that the diaphanous condition efj, = 
constant is an important simplifying element in the anal- 
ysis. If this condition were given up, the problem would 
be very difficult to solve. As mentioned also in jsH], the 
effect is the same as that encountered in the Casimir the- 
ory of a solid ball: the condition of diaphanousness causes 
the divergent terms to vanish [gO^. Analogously, when 
calculating the Casimir energy for a piecewise uniform 
string, the same effect turns up. If the velocity of sound 
(in this case sound replaces light) is the same {— c) in in 
the different pieces of the string, then the theory works 
smoothly [t^I- If this condition is relaxed, the problem 
becomes in practice intractable. 
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leaving us with the radial reduced Green's function equa- 
tion. 



1 d d 2 V 
p op op 



g{p,p')^U{p-p'). (A4) 



The latter, for a Dirichlet arc at p = a, has the familiar 
solution, 



gip.p') = I,Ai^P<)K,Akp>) ~ It,{Kp)I,,{Kp') ^/j-'^^^ , 



p,p <a, 



(A5a) 



9{P:P') = h{K.P<)Ky{Kpy) - K„{Kp)Ky{Kp') 



(A5b) 



p,p > a. 



The azimuthal eigenvalue v is determined by Eq. (jA3p . 
For the wedge 5-function potential (jAip it is easy to de- 
termine V by writing the solutions to Eq. (|A3[1 as linear 
combinations of e^'"^^, with different coefficients in the 
sectors \9\ < a/2 and tt > \d\ > a/2. Continuity of the 
function, and discontinuity of its derivative, are imposed 
at the wedge boundaries. The four simultaneous linear 
homogeneous equations have a solution only if the secular 
equation is satisfied: 



= D{v) = sin^ iy{a - vr) 



1 - 



Sm TTV 



1" T~ I sm 2'KV. 

Ai A2 , 



(A6) 



Because we recognize that the reflection coefficient for a 
single 5-function interface is = (1 + 2w/Ai)~^, so 



APPENDIX A: SEMITRANSPARENT WEDGE 

In this appendix we sketch another way of deriving 
the azimuthal dependence, based on an analogous scalar 
model, in which the wedge is described by a (5-function 
potential, 

V{e) = \i6{e-a/2) + \2&{e + a/2). (Al) 

This has the diaphanous property of preserving the speed 
of light both within and outside the wedge. We can solve 
this cylindrical problem in terms of the two-dimensional 
Green's function G, which satisfies 



ld_ d_ 

p dp^ dp 



v{e) 



G{p, 



p 



-5{p^p')5{e~e'). 



■.p'.o') 

(A2) 



This separates into two equations, one for the angular 
eigenfunction Qu{()) 



92 



(A3) 



^er^^r:^^ = 1 - 



— -, ^1^2^ = — + —, (A7) 

A1A2 Al A2 



we see that this dispersion relation coincides with that in 
Eq. (|1.20p when the reflection coefficient is purely real. 
Note that the v — root of Eq. (|A6|) is spurious and 
must be excluded; there are no = modes for the 
semitransparent wedge. 

Now the full Green's function can be constructed as 

G{x,x') = / ^e-^"(*-*') / *e'fe(--') 
./ 2tt I 2iT 



P^e.(0)e:(0').9.(p,p'),(A8) 



27r 



from which the Casimir energy per length can be com- 
puted from 



g = 



where we have recognized that because the eigenvalue 
equation for v is a, Sturm-Liuoville problem, the integra- 
tion over the eigenfunctions is 27r. As above, we can 
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enforce the eigenvalue condition by the argument princi- noting that ^ In D{iri) is odd, which then yields the ex- 
pie, so we have the expression after again converting to pression 
polar coordinates as in Eq. (|2.10p . 



1 



an K "-V \ ~r' 

J-oo \dri 



In D{iri) 



dppgtriip, p)- 



dx : 



drj — IxiDiirf) j arctan 



\drj 



dpP9zr,{p, P) 



a d 
2k, dua 



lii[Ii^{Ka)KirjiKa)], (All) 



(A13) 



(AlO) where 
Further, we must subtract off the free radial Green's func- 
tion without the arc at r = a, which then implies 



2 sin up, 



[I^.{x)-I^^{x)], (A14a) 



2 sinTT/i 



[lM+I-^{x)], (A14b) 



as well as remove the term present without the wedge where both Lin{x) and Kij^{x) are real for real -q and x, 
potential: and 



A1A2 D{v) 



(A12) 



sinh?77r 

Iir,[x) ^ [Liy^{x) - iKi,^{x)\. (A15) 



leaving us with an expression for the Casimir energy anal- 



Details of the calculation of the Casimir energy for a 



ogous to Eq. (|2.12p . This can be further simplified by semitransparent wedge will appear elsewhere. 
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